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1.  Abstract 

Over  the  last  five  years  we  have  fulfilled  our  contributions  to 
analytical,  numerical  and  experimental  aspects  of  oblique-shock-wave 
reflections  in  pure  and  dusty  gases.  One  of  the  major  contributions  has 
been  the  validation  of  numerical  simulations  with  interferometric  data 
as  long  as  strong  viscous  interactions  were  not  present.  This  has  lent 
confidence  to  the  numerical  simulations  of  blast-wave  phenomena  generated 
by  chemical  or  nuclear  explosives,  where  experimental  data  is  not  avail¬ 
able.  It  has  justified  the  expenditure  of  large  sums  of  money  to  do 
such  numerical  simulations.  Our  work  has  also  pointed  to  using  the 
Navier-Stokes  equations,  rather  than  the  present  Euler  equations,  for 
numerical  simulations  with  significant  viscous  effects  and  to  compare  with 
our  existing  interferometric  data.  The  best  simulations  to  date  have 
been  achieved  by  Dr.  H.  Glaz  and  Dr.  P.  Col  el  la  (see  subsequent  references). 

Another  important  contribution  has  been  the  prediction  and  verification 
of  the  existence  of  Terminal  Double  Mach  Reflection  (TDMR).  We  will 
have  more  to  say  about  this  subject  after  Mr.  J.  Urbanowicz  completes 
his  M.A.Sc.  thesis,  in  the  very  near  future,  and  at  the  7th  Mach 
Reflection  Symposium,  scheduled  for  9-11  June,  1987,  at  Albuquerque, 
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N.M.,  under  the  Chair  of  Charles  Needham  of  S  . 

We  have  published  six  papers  by  Miura  and  Glass  on  dusty-gas  flows 
in  the  Proceedings  of  the  Royal  Society  (London)  and  a  seventh  has  been 
submitted  as  well,  on  many  aspects  of  these  flows.  They  have  laid  the 
analytical  foundations  for  investigating  these  types  of  flows.  In  addition, 
we  have  solved,  for  the  first  time,  the  compressible  laminar  boundary 
layer  flows  for  a  dusty  air  mixture  for  a  semi-infinite  flat  plate  and 
for  the  dusty  gas  boundary  layer  induced  by  a  moving  shock  wave  on  the 


walls  of  a  shock  tube. 


Unfortunately,  our  dusty-gas  flow  experiments  have  not  kept  pace 
with  the  analytical  work.  Although  we  have  a  very  fine  shock-tube  facility 
for  this  purpose,  which  has  been  run  as  a  classical  shock  tube  for 
pressure-gauge  calibrations,  we  have  not  succeeded,  as  yet,  in  injecting 
dust  to  produce  a  homogeneous  air-dust  mixture/  We  have  tried  several 
unsuccessful  techniques  and  consulted  a  number  of  laboratories  on  this 
problem,  it  looks  as  we  will  have  to  come  up  with  a  unique  solution 
for  our  facility.  We  are  hopeful  of  doing  this  in  the  near  future. 

Finally,  our  contributions  to  oblique-shock-wave  reflections  will 
be  summarized  in  an  abbreviated  version  of  my  AIAA  distinguished  Dryden 
Lecture  in  Research  in  the  AIAA  Jl.  entitled,  "Some  Aspects  of  Shock- 
Wave  Research",  which  is  due  to  appear  soon. 

2.  Details  of  Accomplishments  and  Publications  in  Journals  and  UTIAS 
Reviews,  Reports  and  Technical  Notes 

a)  List-  of  Publications  from  1982  to  1987  in  Journals,  Over  the  Past  Five  Years 
(additional  publications  are  forthcoming  in  1957) 

1.  Glass,  I.  I.,  "Some  Aspects  of  Shock  Wave  Research",  AIAA  J.,25,2,  pp. 
214-219,  (1987). 

2.  Hu,  T.  C.  J.  and  Glass,  I.  I.,  "Pseudostationary  Oblique  Shock-Wave 
Reflections  in  Sulphur  Hexafluoride  (SFg):  Interferometric  and  Numerical 

3.  Miura,  H.  and  Glass,  I.  I.,  "Oblique  Shock  Waves  in  a  Dusty-Gas  Flow 
Over  a  Wedge",  Proc.  Soc.  Lond.  A  408,  pp.  61-78  (1986). 

4.  Hu,  T.  C.  J.  and  Glass,  I.  I.,  "Blast  Wave  Reflection  Trajectories  from  a 
Height  of  Burst",  AIAA  J.,  Vol .  24,  No.  4,  pp.  607-610,  (1986). 

5.  Glaz,  H.  M.,  Glass,  I.  I.,  Li,  J.-C.  and  Walter,  P.  A.,  "Interaction  of 
Oblique  Shock-Wave  Reflections  in  Air  and  CO2  with  Downstream  Obstacles", 
in  Shock  Waves  and  Shock  Tubes,  edited  by  D.  Bershader  and  R.  Hanson, 
Proceedings  Fifteenth  International  Symposium  on  Shock  Waves  and  Shock 
Tubes,  Stanford  University  Press.  Stanford,  CA,  (1986). 

6.  Hu,  T.  C.  J.  and  Glass,  I.  I.,  "An  Interferometric  and  Numerical  Study 

of  Pseudo-Stationary  Oblique-Shock-Wave  Reflections  in  Sulfur  Hexafluoride 
(SFc)",  Fifteenth  International  Symposium  on  Shock  Waves  and  Shock  Tubes, 
(1986). 
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7.  Li,  J.-C.  and  Glass,  I.  I.,  "Collision  of  Mach  Reflections  with  a  90-Degree 
Ramp" ,  Fifteenth  International  Symposium  on  Shock  Waves  and  Tubes,  (1986). 

8.  Glaz,  H.  M. ,  Walter,  P.  A.,  Glass,  I.  I.  and  Hu,  T.  C.  J.,  "Oblique  Shock 
Wave  Reflections  in  SF6:  A  Comparison  of  Calculation  and  Experiment", 

in  Dynamics  of  Explosions,  Proceedings,  Tenth  International  Colloquium 
on  Dynamics  of  Explosions  and  Reactive  Systems,  Berkeley,  CA,  edited  by 
J.  R.  Bowen,  J.-C.  Leyer  and  R.  I.  Soloukhin,  published  by  AIAA  Inc., 
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9.  Shirouzu,  M.  and  Glass,  I.  I.,  "Evaluation  of  Assumptions  and  Criteria  in 
Pseudostationary  Oblique  Shock-Wave  Reflections",  Proc.  R.  Soc.  Lond. 
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10.  Kleiman,  J.,  Chan,  P.  W.  C.,  Salansky,  N.  and  Glass,  I.  I.,  "Chemical- 
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J.  Appl .  Phys.  59(6),  pp.  1956-1961,  (1986). 

11.  Miura,  H.,  Saito,  T.  and  Glass,  I.  I.,  "Shock-Wave  Reflection  from  a  Rigid 
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Conclusions 
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continue  in  the  future  from  the  three  respective  aoencies. 

We  in  turn  have  tried  our  best  to  provide  good  value  .'or  the  dollar  received 
in  the  form  of  recognized  research  contributions  in  our  UT I  AS  Reports,  Technical 
Notes  and  Reviews,  as  well  as  in  the  open  literature  puDlications  in  the  Proc. 

R.  Soc.  (London),  J.  Fluid  Mech.,  Phys.  Fluids,  AIAA),  .1.  Appl .  Phys.  and 
the  various  Mach-Refl ection  and  Shock-Tube  Proceedings. 
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Nomenclature 
a  =  speed  of  sound 

C  =  wedge  corner 

CMR  =  complex-Mach  reflection 

DMR  =  double-Mach  reflection 

e  =  specific  internal  energy 

E  =  total  specific  energy,  '/i(u2  +  v2)  +  e 

h  -  flow  enthalpy 

H  =  pseudostationary  total  enthalpy, 

=  + 1?)  +  h 

/  =  incident  shock  wave  in  shock  tube 

K  =  kink  in  CMR  reflected  wave  R 

M.M'  =  first  and  second  Mach  stems  in  DMR 

ffi  =  self-similar  Mach  number  \M),  =  (u2  +  P)/a2 
MR  =  Mach  reflection  (SMR,  CMR,  DMR,  TDMR) 

M,  =  incident  shock  wave  (/)  Mach  number 
p  =  static  pressure 

P  =  reflection  point;  point  where  RR -CMR  and 

CRM  -  DMR  lines  meet 
R  =  gas  constant 

R, R’  =  first  and  second  reflected  shock  waves  in  DMR 

RR  =  regular  reflection 

S, S'  =  first  and  second  slipstreams  in  DMR 

SMR  =  single-Mach  reflection 

/  =  time 

T  *  temperature 

TDMR  =  terminal  double-Mach  reflection,  where  x'  =  0 

T,  T'  =  first  and  second  triple  points  in  DMR 

u  =  (u,v)  velocity  field 

u  =m-{  pseudostationary  velocity  component  in  the 

x  direction 

v  =v-*i  pseudostationary  velocity  component  in  the 

y  direction 

u,  v  =  velocity  components  in  the  x  and  y  directions 
x.y  =  Cartesian  coordinates 

a  =  velocity-deflection  angle  in  slipstream  shear  layer 

y  =  specific  heat  ratio 

6  =  angle  between  incident  (/)  and  ref  cted  (R) 

shock  waves 

5*  =  boundary-layer  displacement  thickness 

t  =  Mach-stem  curvature  angle 

i,  =O-y0)/(/-r0)  pseudostationary  coordinate 


8d  =  boundary-layer  displacement  effect  on  the  dif¬ 

ference  in  the  flow  angles  of  the  incident  and 
reflected  shock  waves 
02*,  =  flow  detachment  angle 

0.  =  wedge  angle 

6‘w  =  effective  wedge  angle,  =  0»  +  x 

{  =  pseudostationary  coordinate,  =  (x  -  x0)/ (/  - 10) 

p  =  flow  density 

4  =  wave  angle 

X,  x '  =  first  and  second  triple-point  trajectory  angles 


Introduction 

HOCK-WAVE  phenomena  on  Earth1-1  are  common  and 
important  occurrences  in  our  lives.  For  example,  we  hear 
them  from  the  cradle  to  the  grave  as  terrorizingly  crashing  or 
rumbling  thunder  following  lightning  discharges.4-7  It  is 
estimated  that  1000  thunderstorms  occur  on  Earth  at  any  mo¬ 
ment.  In  the  primordial  state  of  our  Earth,  such  shock  waves 
are  credited  as  being  enormously  more  efficient  in  creating  the 
building  and  replicating  blocks  of  life,  such  as  ATP,  RTP,  and 
DNA,2'*  than  the  radiation  from  the  sun.  Yet,  paradoxically 
enough,  all  life  is  threatened  today  by  the  radiation  and  the 
very  shock  waves  that  might  be  unleashed  by  a  thermonuclear 
war.M1 

The  lethal  overpressures  and  winds  induced  by  shock  waves 
are  now  well  understood,1211  even  though  our  eyes  cannot  see 
the  shock  waves  generated  by  lightning,  firearms,  or  explo¬ 
sions  in  air,  since  air  is  a  transparent  gas.  Special  optical 
methods  (interferometry,  schlieren,  and  shadow  photography) 
must  be  used  to  make  shock  waves  visible.14  The  more  intense 
the  explosion  or  rapidity  and  size  of  energy  release,  the  higher 
is  the  velocity  of  the  resulting  shock  wave,  and  its  attendant  in¬ 
creases  in  pressure,  density,  and  flow  velocity,  and  the  greater 
is  the  possible  devastation.  The  shock  wave  becomes  less 
destructive  with  time  or  distance  as  it  compresses  and  heats  the 
enveloped  air,  thereby  dissipating  its  energy,  until  it  becomes  a 
harmless  sound  wave.1214 

At  sea  level,  where  the  mean  free  path  is  small  (6.6  x  10* 6 
cm),  the  thickness  of  the  shock  front  is  about  tenfold  this 
value.  Consequently,  it  is  quite  incredible  how  rapidly  the 
shock  wave  changes  from  its  quiescent  ambient  conditions  in 
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front,  to  one  of  destruction  behind  it  through  molecular  colli¬ 
sions.  In  outer  space,  where  the  mean  free  path  is  very  targe, 
the  shock  thickness  can  be  thousands  of  kilometers,  and  the 
various  changes  through  it  are  much  more  gradual. 

Aside  from  the  lightning-thunder  process,  shock  waves  on 
Earth  are  also  generated  by  intense  volcanic  eruptions  and 
meteor  impact.  The  recent  eruption  of  Washington’s  Mt.  St. 
Helens  is  estimated  to  have  had  an  energy  release  equivalent  to 
50  Mtons  of  TNT  or  about  3000  bombs  of  the  type  used  to 
destroy  Hiroshima.  The  eruption  of  the  Krakatoa  volcano  in 
1883  in  Indonesia  is  estimated  to  have  had  an  energy  release  of 
5000  Mtons  of  TNT  or  100  times  as  powerful  as  Mt.  St. 
Helens.  Perhaps  the  greatest  eruption  in  recorded  history  oc¬ 
curred  in  Tambora,  Indonesia,  in  1815,  with  an  energy  release 
of  about  tenfold  that  of  Krakatoa.  The  airborne  ash  caused  “a 
year  without  a  summer”  in  1816  and  12,000  people  were 
Idled.15 

If  a  meteor  of  the  size  that  produced  the  Arizona  Crater 
ever  impacted  in  a  big  city,  it  could  be  as  devastating  as  a 
multimegaton  bomb.16  An  energy  release  equivalent  to  15 
Mtons  TNT  is  estimated  to  have  been  released  on  impact  over 
50,000  years  ago.  Yet,  there  is  geological  evidence  of  craters 
formed  millions  of  years  ago  from  ap  energy  release  of 
50  million  Mtons  of  TNT,  caused  by  asteroid  impact.  Such 
asteroids  orbit  the  sun  between  Mars  and  Jupiter.  A  few  orbit 
the  sun  to  pass  inside  the  orbit  of  the  Earth.  However,  not  all 
are  known.  It  is  estimated  that  perhaps  800  ±  300  such  Apollo 
asteroids,  larger  than  1  km,  may  remain  undetected.  An  im¬ 
pact  would  release  an  energy  of  about  100,000  Mtons  of  TNT 
and  leave  a  crater  20  km  wide.  As  a  matter  of  fact,  a  news  item 
[Washington  (Reuter)  14.2.81  ]  reports  that  NASA  urged  the 
U.S.  government  to  organize  “Project  Spacewatch”  to  keep 
track  of  such  asteroids  and  meteors  and,  if  necessary,  to  send 
spacecraft  armed  with  hydrogen  bombs  to  deflect  those  that 
may  menace  to  impact  the  Earth. 

Yet,  we  need  not  associate  shock  waves  only  with  destruc¬ 
tion.  Our  modern  industrial  society  could  not  have  been  built 
without  the  use  of  explosives  that  produce  shock  waves.  Ex¬ 
plosives  are  used  in  building  highways,  canals,  tunnels,  har¬ 
bors,  railroads,  and  subways;  in  mining,  agriculture,  engineer¬ 
ing;  and  in  the  space  programs  for  precise  timing,  activating, 
and  cutting  requirements.1  Recently,  shock  waves  produced 
by  a  high-voltage  discharge  or  5-10  mg  of  lead  azide  (PbNt) 
have  been  applied  successfully  to  medicine  in  the  breaking  up 
of  large  kidney  stones  and  bladder  stones  to  relieve  suffering 
and  possibly  death.1711  More  recently,  lasers  have  been  used 
as  an  energy  source  for  the  same  purpose  with  some  promising 
results.19  It  is  estimated  that  50,000  successful  operations  have 
been  performed  using  the  spark-discharge  machine  described 
in  Ref.  17,  which  is  quite  expensive.  The  cost  factor  (several 
million  dollars)  is  important  for  use  in  the  Third  World,  where 
the  occurrence  of  such  stones  are  more  common.  It  partly 
motivated  the  work  in  Refs.  18  and  19  to  reduce  costs  by  a  fac¬ 
tor  of  10,  at  least.  Basically,  the  three  systems  are  similar  in 
that  the  energy  discharge  (spark,  explosive,  laser)  takes  place 
at  the  focus  of  an  ellipsoidal  mirror  and  is  refocused  at  the  sec¬ 
ond  focus  located  on  the  kidney  stone,  as  a  powerful  implo¬ 
sion  wave,  which  breaks  up  the  obstructing  stone.  It  may  take 
many  shots  to  break  up  a  stone.  In  the  case  of  bladder  stones, 
the  explosive  is  used  directly  against  the  stone. 

In  space,  the  best  example  of  shock-wave  phenomena  is  the 
magnetosphere1  flow  about  planets  in  our  solar  system.  Our 
Earth  is  an  excellent  example.  It  becomes  a  spherical -flow 
model  in  a  solar  wind  at  hypersonic  speeds  (M ~  13).  The  solar 
wind  is  composed  of  ionized  hydrogen,  helium,  some  other 
heavier  nucleons  and  electrons.  This  plasma  is  electrically 
neutral.  It  blows  with  a  supersonic  velocity  of  400  km/s  with  a 
bulk  number  density  of  about  10  particles/cm’  having  an 
average  kinetic  temperature  of  about  105  K.  The  interaction  of 
the  ionized  particles  with  the  Earth’s  Jipole  magnetic  field 
distorts  it  into  a  streamlined  shape  w  jth  the  magnetopause  as  a 
boundary,  that  is,  where  the  ionized  gas  can  no  longer 


penetrate  the  Earth’s  magnetic  field.  The  collisionless  bow 
shock  in  from  of  the  line  of  symmetry  is  about  30  Earth  radii 
from  the  Earth's  center.  The  stagnation  point  on  the 
magnetopause  is  about  20  Earth  radii  from  the  Earth’s  center 
on  this  line.  The  standoff  distance  of  the  bow  shock  from  the 
stagnation  point  is  about  10  Earth  radii  (about  64,000  km). 
The  bow  shock  itself  is  many  thousands  of  kilometers  thick. 
The  subsonic  plasma  behind  the  bow-  shock  probably  has  a 
tenfold  increase  in  temperature  to  about  106  K.  Near  the  Earth 
the  magnetopause  cavity  contains  the  famous  Van  Allen  radia¬ 
tion  belts.  The  magnetopause  cavity  resembles  a  comet,  and 
the  tail  may  extend  several  million  kilometers  downwind. 

Similar  flows  can  exist  around  a  star.  A  recent  relevant 
study  of  this  phenomenon  is  given  in  Ref.  20.  Of  the  many 
shock-wave  phenomena  in  space,  one  of  the  most  interesting  is 
the  existence  of  shock  waves  in  the  spiral  arms  of  galaxies21 
such  as  our  own  Milky  Way.  Densities  can  increase  by  about 
an  order  of  magnitude,  making  possible  the  birth  of  young 
blue  stars. 

Pseudostationary  Oblique  Shock-Wave  Reflections 
.  This  problem  was  first  investigated  by  E.  Mach  in  1878, 22 
and  has  since  occupied  many  scientists  and  engineers,  and 
where  my  students  and  I  have  made  some  worthwhile  con¬ 
tributions  over  the  past  decade.  The  problem  is  known  as 
“pseudostationary  oblique  shock-wave  reflections.”  Al¬ 
though  much  more  is  known  about  this  problem  now  than 
when  Mach  started  it  all,  there  is  still  a  great  deal  to  learn.  It  is 
of  interest  for  two  main  reasons:  first,  the  problem  involves  a 
great  deal  of  gasdynamics  of  perfect  and  imperfect  flows 
(frozen  and  equilibrium,  inviscid  and  viscous);  second,  our  ex¬ 
perimental  (interferometric)  shock-tube  data  have  provided 
computational  fluid  dynamicists,  who  are  working  on  explo¬ 
sion  dynamics,  with  the  means  of  checking  their  numerical 
simulations  for  accuracy.  The  latter  data  provide  engineers 
with  pressure,  velocity,  and  other  flow  properties  required  to 
build  explosion-resistant  structures.  To  date,  such  simulations 
have  proved  to  be  extremely  good,  considering  that  the 
numerical  solutions  solve  only  the  inviscid  Euler  equations, 
whereas  the  experimental  data  come  from  flows  that  are 
viscous  and  can  have  real-gas  effects  as  well. 

Analysis 

Unlike  linear  waves  encountered  in  acoustics  and  light, 
shock  waves  are  nonlinear.  For  example,  the  reflected  wave 
angle  is  different  from  the  angle  of  incidence,  and  the  physical 
quantities  across  shock  waves  such  as  pressure,  density, 
temperature,  and  flow  velocity  are  highly  nonlinear.  The 
simplest  method  of  studying  pseudostationary  oblique  shock- 
wave  reflections  is  in  a  shock  tube.  Here  a  planar  shock  wave 
M,  collides  with  a  sharp  compressive  wedge  of  fixed  angle  0„, 
and  gives  rise  to  regular  reflection  (RR)  and  three  types  of 
Mach  reflections  (MR)  —  namely,  single-  (SMR),  complex- 
(CMR),  and  double-Mach  reflection  (DMR)  —  as  shown  in 
Figs.  1  and  2.  (The  various  symbols  arc  defined  on  the 
figures.)  It  is  seen  that  RR  consists  of  two  shock  waves 
meeting  at  the  reflection  point  P.  However,  MR  consists  of 
three  shock  waves  meeting  at  the  triple  point  T.  It  has  been 
shown2’  25  that  such  flows  are  self-similar,  i.e. ,  they  look  the 
same  with  time  except  for  an  increased  scale.  Consequently, 
the  Mach  stem  grows  with  time  (or  distance)  along  the  triple¬ 
point  trajectory  angle  x-  It  may  also  be  mentioned  that  this 
angle  represents  an  effective  wall  for  the  incident  shock  wave  / 
to  reflect  from  and  produce  a  reflected  shock  wave  R.  This 
provides  the  subsequent  use  of  the  effective  wedge  angle 
61  =  8,  +  x  as  a  useful  parameter. 

It  is  seen  that  regular  reflection  consists  of  an  incident  shock 
wave  /  and  a  reflected  shock  wave  R,  which  may  be  attached  at 
the  wedge  corner  C,  or  detached  from  it.  The  standoff 
distance  j,  divided  by  the  distance  L,  which  the  incident  shock 
wave  /  has  traveled  from  the  corner,  remains  constant.  The 


(d)  x’>  0,  w'<  0 


Ftf.  2  Typo  of  double-M«ch  reflections  in  pseudostationary  flow,  a)  T’ 
lies  iboul  x<  b)  T'  lies  along  *,  c)  T'  li«  below  x,  «•)  terminal  double-Mach 
reflection  (TDMR)  where  T  ’  lies  along  wedge  surface,  M'  has  vanished,  and 
R'  Is  a  regular  reflection. 


reflected  shock  wave  is  curved  if  the  flow  behind  it  is  subsonic, 
and  it  may  have  a  straight  portion  near  the  reflection  point  if 
the  flow  there  is  supersonic  and  uniform  (see  Figs,  la  and  1c). 

Mach  reflection  (Figs.  I b- Id),  on  the  other  hand,  consists 
of  incident  /  and  reflected  R  shock  waves,  a  Mach  stem  M, 
and  a  slipstream  S.  The  three  shocks  meet  at  the  triple  point  T, 
and  the  slipstream  S  separates  states  (2)  and  (3),  even  though 
they  have  the  same  pressure  and  flow  inclination.  SMR  has  a 
curved  reflected  wave.  CMR  has  a  kink  in  the  reflected  wave, 
which  is  caused  by  an  interacting  compression  wave.  In  DMR, 
the  kink  becomes  sharp  and  a  second  Mach  shock  M'  is 
formed  that  finally  extends  to  the  slipstream  S.  As  a  result,  a 
second  triple  point  T'  and  a  second  slipstream  S'  are  also 
formed.  In  this  case,  the  reflected  portion  R  is  straight  and  the 


flow  behind  it  is  supersonic.  The  reflected  shock  R  in  SMR 
and  R'  in  all  other  Mach-reflection  (MR)  cases  can  also  be  at¬ 
tached  or  detached  (see  Figs,  lb- Id  and  2b-2d). 

The  question  arises,  given  the  initial  conditions  of  the  gas 
p0,  T0,  y,  the  shock  Mach  number  M,,  and  the  wedge  angle 
6W,  can  one  predict  the  resulting  reflection  a  priori?  The 
answer  is  “yes."  This  can  be  obtained  from  the  following 
analysis. 

It  has  been  noted  previously  that  such  two-dimensional 
wedge  flows  are  self-simitar.  Consequently,  instead  of  three 
independent  variables  x,  y,  and  r,  the  reflection  phenomenon 
can  be  described  in  terms  of  two  independent  variables  x/r  and 
y/t.  The  reflection  point  P,  in  regular  reflection,  and  the  triple 
point  T,  in  Mach  reflection,  are  chosen  as  specific  points  mov- 
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ing  at  constant  velocity  with  respect  to  the  wedge  corner  C.  By 
attaching  frames  of  reference  to  these  points,  the  reflections 
become  pseudosteady  (Fig.  I).  Consequently,  the  steady-flow 
equations  of  motion  can  be  applied  to  each  shock  wave  in 
turn,  such  that 

Continuity: 

P/u<  =P/U/  sin(*,  -0y )  (1) 

Tangential  momentum: 

p,  tan*,  tan(*,-0y)  (2) 

Normal  momentum: 

P,  +  p,uj  sin2*,  =Pj  +  Pju]  sin2(*,  -0, )  (3) 

Energy: 

h,  +  Viu]  sin2*,  =  hj  +  Vi u)  sin2(*, -0y )  (4) 

where  p,  is  the  flow  density  of  the  initial  state  in  front  of  the 
shock  wave  and  p,  is  the  flow  density  in  the  final  state  behind 
the  shock  wave,  the  same  for  the  flow  velocity  u,  pressure  p, 
specific  flow  enthalpy  h,  wave  angle  *,  and  flow-deflection 
angle  0.  An  equation  of  state 

p  =  PRT  (5) 

is  also  used  for  a  perfect  gas  or  p=p{p,T)  for  an  imperfect 
gas.  In  addition,  the  boundary  condition  for  regular  reflection 
must  be  used  so  that  the  flow  is  parallel  to  the  wedge  surface 
(Fig.  1) 

0i  =«2  (6) 

For  Mach  reflection,  across  the  contact  surface  or 
slipstream,  the  pressure  must  be  continuous 

Pi=P>  C?) 

and  the  flow  deflections  must  also  be  continuous  or 

0j  =  0i  -  63  (8) 

It  has  been  shown  by  von  Neumann26-27  and  verified  by 
many  experimentalists  that  in  pseudostationary  flow  a  regular 
reflection  will  undergo  transition  to  a  Mach  reflection  when 
02  =  02».  that  is,  when  the  wedge  angle  0„  produces  the  max¬ 


imum  possible  flow  deflection  0^.  This  angle  is  also  known  as 
the  flow-detachment  angle.  For  further  details,  see  Refs.  14 
and  26-28. 

As  noted  previously,  Mach  reflection  (Fig.  1)  can,  in  turn, 
be  divided  into  SMR,  CMR,  and  DMR;  and  DMR  can  be  sub¬ 
divided  into  four  subtypes  (Fig.  2),  depending  on  whether  or 
not  the  second  triple-point  trajectory  angle  * '  is  above,  equal 
to,  or  below  the  first  triple-point  trajectory  angle  x-  Finally, 
there  is  the  terminal  double-Mach  reflection  (TDMR)  case, 
where  x '  =  0.  That  is,  the  second  triple  point  lies  on  the  wedge 
surface.  The  second  Mach  shock  disappears  and  the  reflected 
shock  wave  reflects  as  a  regular  reflection.  Such  reflections  oc¬ 
cur  readily  in  gases  with  low  values  of  y,  such  as  Freon- 1225  or 
SFt.M  They  are  not  possible  in  perfect  monatomic,  diatomic, 
or  triatomic  gases,  although  they  are  possible  when  these  gases 
are  imperfect  and  have  a  low  effective  y  (see  Fig.  2d). 

It  is  assumed  that  SMR  becomes  a  CMR  when  the  flow  with 
respect  to  the  first  triple  point  T  becomes  sonic,  i.e.,  M2r=  I. 
This  assumption  is  based  on  experimental  observation  and 
computational  simulation.2*' 25  The  same  heuristic  approach  is 
taken  for  CMR -DMR  transition,  and  the  condition  is  ap¬ 
plied  to  the  Vink  or  second  triple  point  K  or  T  ’  or  =  1 .  A 
change  of  curvature  occurs  in  the  reflected  wave  /?,  which 
steepens  into  a  sharp  second  triple  point  T'.  This  marks  the 
formation  of  the  second  Mach  shock  M',  and  a  straight 
reflected  shock  R,  before  M‘,  and  a  curved  reflected  shock 
R‘ ,  after  M' .  The  flow  behind  R  is  supersonic.  It  should  be 
noted  that  M’  grows  in  strength  and  length  until  it  comes  close 
to  the  slipstream  5.  Ideally,  it  should  terminate  there.  Then  it 
raises  the  problem:  How  can  the  contact  front  remain  stable 
with  a  sharp  pressure  gradient  applied  to  it?  Interferometri- 
cally,  the  Mach  shock  appears  to  go  around  the  rolled-up 
slipstream;  in  this  case,  as  a  compression  wave.36-5*  This  is 
also  verified  by  numerical  simulation  of  the  interferograms  in 
question  for  DMR  in  different  gases. J3  M 

If  the  preceding  equations  and  conditions  are  applied  and 
solved,  it  is  possible  to  obtain  a  plot  in  the  (M,,  0W,  or  0^) 
plane,2*  delineating  the  various  regions  of  RR,  SMR,  CMR, 
and  DMR  and  their  transition  lines,  as  shown  in  Figs.  3  and  4, 
for  example,  for  perfect  and  imperfect  air  in  vibrational, 
dissociationa),  and  ionizational  equilibrium,  over  a  very  large 
initial  shock  Mach  number  range  l<A/,<20.  It  can  be  seen 
that,  for  Ms>  10,  the  curves  do  not  change  for  perfect  air,60 
but  undulate  for  equilibrium  air  as  new  degrees  of  freedom  arc 
excited. 

By  assuming  that  the  Mach  stem  is  straight  and  perpen¬ 
dicular  to  the  wedge  surface,  from  the  triple  point  T,  it  is 
possible  to  obtain  an  expression  for  the  triple-point  trajectory 
angle  x.41  as  well  as  the  second  triple-point  trajectory  angle 
x'.42  Unfortunately,  the  Mach  stem  is  usually  curved.  It  is 
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Fif(.  Sa  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  (A/,,#,)  plane  for  frozen  argon  (>  =  5/3)  and  ex- 
perimrntal  data.  For  SMR -CMR  transition  boundary,  solid  line 
represents  M 2r  =  I  and  dashed  line  represents  4  =  90  deg.  Experimen¬ 
tal  data:  oRR.  oSMR,  A  (  MR,  •  DMR  (Ref.  31). 


60  r  A  ♦  **  -<* 

RR 


Fig.  Sb  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  <A/,,#»)  plane  for  frozen  argon  <y  =  5/3)  and  ex¬ 
perimental  data.  For  SMR -CMR  transition  boundary,  solid  line 
represents  M2r  =  1  and  dashed  line  represents  4-90  deg.  Experimen¬ 
tal  data:  ORR.  OSMR,  a  (MR,  a  l)MR  (Ref.  31). 


concave  at  the  lower  shock  Mach  number  A/,,  for  significant 
stem  heights,  and  convex*’  at  higher  Mach  numbers.  This  in¬ 
troduces  an  error  in  determining  \  and  x ' .  especially  when  x  is 
small.45 

Comparison  with  Experimental  Results  in  the  M,,  f,,  or  t‘w  Plane 

The  interferometric,  schlieren,  and  shadowgram  data  from 
the  Institute  for  Aerospace  Studies  were  obtained  in  the 
10  x  18-cm  Hypervelocity  Shock  Tube.  The  major  instrument 
for  this  purpose  was  a  23-cm-diam  field-of-view  Mach- 
Zehnder  interferometer  with  a  giant-pulse  Q-switched  ruby 
laser  with  a  frequency  doubler  at  6943  and  3471.5  A,  having 
an  exposure  time  of  15  ns.  The  latter  was  short  enough  to 
photograph  the  strongest  shock  wave  without  blurring. 

Figures  5-9  show  comparisons  of  the  analytical  domains 
and  transition  boundaries  for  RR,  SMR,  CMR,  and  DMR 
with  experiments  in  argon,5'  oxygen,*1  nitrogen,50  air,5*  car¬ 
bon  dioxide,51  and  sulfurhexafluoride.51  In  addition,  the 
detailed  results  obtained  by  Ikui  et  al.2!  for  air,  C02,  and 
Freon- 1 2  are  shown  and  discussed. 

An  examination  of  Fig.  5a  shows  the  experimental  results 
for  RR  and  MR  in  Ar  on  an  (A/,,0„)  plane.  In  our  early  ex¬ 
periments,  the  major  effort  was  to  show  that  the  various 
regimes  were  predictable  from  the  foregoing  analysis.  As  more 
experience  was  gathered,  it  become  clear  that  the  various 
transition  lines  were  the  sensitive  lines  to  be  tested.  Conse¬ 
quently,  the  various  analytical  regimes  look  reasonably  well 
validated.  Three  DMR  runs  (1  at  50  deg  and  2  at  40  deg)  lie  in 
the  CMR  regime.  All  of  the  other  experiments  are  in  their  ap¬ 
propriate  regimes.  It  has  been  found*5  that,  along  with 
A/27-=I  for  SMR-CMR,  another  condition  that  must  be 
satisfied  is  6>90  deg.  On  all  graphs,  the  dashed  line  represents 
6  =  90  deg  and  its  lower  branch  lies  above  the  solid  line  for 
SMR-CMR,  where  M:r  =  1,  The  upper  branch  of  the  dashed 
line  lies  below  the  shock  line,  where  M:t#  1,  yet.  It  can  be 
seen  that  two  points  at  0„  =  40  and  20  deg  are  improved  to  lie 
in  CMR. 

When  these  results  are  plotted  on  an  plane  (Fig. 

5b),  which  is  more  accurate  since  a  relation  for  x  is  not  re¬ 
quired,  paradoxically,  a  poorer  agreement  is  obtained  for  Ar. 
Now  four  DMR  runs,  one  at  6W  =  50  deg  and  three  at  6,  -  45 
deg,  definitely  lie  in  the  CMR  regime  and  two  other  CMR’s 
are  borderline  runs;  one  CMR  run  lies  in  the  SMR  regime  at 
A/,-4,  all  other  runs  are  in  their  appropriate  regimes.  Since 
no  runs  were  made  along  the  RR  — MR  transition  line,  little 
can  be  said  about  the  "von  Neumann  paradox"  (where  RR 
persists  into  the  MR  regime)  or  possible  boundary-layer  ef¬ 
fects.  No  real-gas  effects  occur  in  Ar  for  the  pressures  and 
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Fif .  6a  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  ( M,,tw )  plane  for  frozen  N2  and  02  (y  =  7/5)  and 
experimental  data.  For  SMR -CMR  transition  boundary,  solid  lint 
represents  M22-  =  1  and  dashed  line  represents  6  =  90  deg.  Experimen¬ 
tal  data:  oRR.  oSMR,  a  CMR,  aDMR  (N2— Ref.  30);  a  SMR, 
a  CMR,  o  DMR  (02— Ref.  41). 


Fig.  6b  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  ( M,.tw )  plane  for  frozen  N2  and  02  (y  =  7/5)  and 
experimental  data.  For  SMR -CMR  transition  boundary,  solid  line 
represents  M2r  =  1  and  dashed  line  represents  4  =  90  deg.  Experimen¬ 
tal  data:  oRR,  oSMR,  a  CMR,  •  DMR  (N2 — Ref.  30);  nSMR, 
a  CMR,  o  DMR  (02— Ref.  41).  NR  is  the  no-reflection  region. 


Fig.  6c  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  ( M,,tw )  plane  for  air  and  experimental  data.  Solid 
lines  are  for  air  as  a  perfect  gas  (y  -  7/5).  For  SMR  -  CMR  transition 
boundary,  solid  line  represents  Mjr  =  1  and  dash-dot  line  represents 
4«90  deg.  Dashed  lines  are  for  Imperfect  air  at  Pt  -  15  Tore  and 
Tt  =  300  K.  Experimental  data:  oRR,  oSMR,  a  CMR,  aDMR 
(Ref.  36). 


Fig.  6d  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  ( )  plane  for  air  and  experimental  data.  Solid 
lines  are  for  perfect  air  (y  =  7/5).  Dashed  lines  are  for  imperfect  air  at 
P,  =  15  Torr  and  Tt  *  300  K.  For  SMR  -CMR  transition  bouodary, 
solid  line  represents  Mir  =  1  and  dashed  line  represents  4  =  90  deg.  Ex¬ 
perimental  data:  oRR,  oSMR,  a  CMR,  •  DMR  (Ref.  36).  NR  is  the 
no-reflection  region. 


Fig.  6e  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  In  the  plane  (or  perfect  air  (y-7/5)  and  ex¬ 

perimental  data.  For  SMR -CMR  transition  boundary,  solid  line 
represents  M1T  » 1  and  dashed  line  represents  4  =  90  deg.  Experimen¬ 
tal  data:  oRR,  oSMR,  a  CMR,  •  DMR  (Ref.  25). 


Mach  numbers  of  these  experiments,  as  the  relaxation  lengths 
are  too  large  and  the  gas  remains  frozen  with  a  y  =  5/3. 

Figures  6a  and  6b  show  the  results  of  N2  and  O-  in  the  (M,, 
8W)  and  (A/,,8*)  planes,  in  Fig.  6a,  the  six  RR  runs  at  the 
RR  — MR  transition  are  only  marginally  below  it.  All  runs  ex¬ 
cept  a  DMR  at  A/,  =  6  lie  in  their  appropriate  regimes.  The  line 
for  4  =  90  deg  is  very  helpful  in  this  case,  because  all  of  the 
SMR  runs  lie  in  their  appropriate  regimes.  The  same  com¬ 
ments  apply  to  Fig.  6b.  It  is  seen  that  as  y  decreases  the  NR 
and  SMR  — CMR  lines  lie  closer  together,  preventing  a  clear- 
cut  distinction  of  the  experiments. 

Figures  6c  and  6d  show  the  results  for  air  on  the  (A/,,8*)  and 
(A/,,8*)  planes.34  It  is  seen  from  Fig.  6c  that  at  the  RR-MR 
transition  line  two  RR  runs  lie  in  the  SMR  regime;  one  RR  run 
lies  in  the  CMR  regime  and  seven  RR  runs  lie  in  the  DMR 
regime.  This  definitely  shows  the  persistence  of  RR  into  the 
MR  regime.  Subsequently,  it  will  be  shown  that  this  is  due  to 
the  shock-induced  boundary  layer  on  the  wedge  surface.  It  is 
also  known  that  the  CMR  — DMR  line  must  curve  and  meet 
the  SMR  — CMR  line  at  point  P  on  the  RR-MR  line,  as  the 
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Fig.  7a  Regions  and  transition  lines  of  regular  (RR)  and  Mack  (MR) 
reflections  in  the  (M,,!.)  plane  (or  CO,  as  a  perfect  gas  <y  ■  1.29) 
and  experimental  data.  For  SMR-CMR  transition  boundary,  solid 
line  represents  Mlr=l  and  dashed  line  represents  J  *  90  deg.  Ex¬ 
perimental  data:  o  RR,  □  SMR,  a  CMR,  a  DMR  (Ref.  37). 


Ms 

Fig.  7b  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  (W,,IJ  plane  for  CO)  as  a  perfect  gas  (>  =  1.29) 
and  experimental  data.  For  SMR-CMR  transition  boundary,  solid 
line  represents  M,r  =  I  and  dashed  line  represents  4  =  90  deg.  Ex¬ 
perimental  data:  oRR,  qSMR,  a  CMR,  o  DMR  (Ref.  37). 
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Fig.  7c  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  In  the  (M,,lw)  plane  for  vibrational  equilibrium  CO,  and 
experimental  data.  For  SMR-CMR  transition  boundary,  solid  line 
represents  Mjf»l  and  dashed  line  represents  4*90  deg.  Dash-dot 
lines  are  for  frozen  CO,  (y  *  7/S).  Experimental  data:  o  RR,  □  SMR, 
a  CMR,  a  DMR  (Ref.  37). 
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Fig.  7d  Regions  and  traasilion  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  In  the  plane  for  vibrational  equilibrium  COt  and 

experimental  data.  For  SMR-CMR  transition  boundary,  solid  line 
represents  M,r*  1  and  dashed  line  represents  4*90  deg.  Dash-dot 
lines  are  for  frozen  CO)  (>  *  7/S).  Experimental  data:  o  RR,  □  SMR, 
a  CMR,  a  DMR  (Ref.  37).  NR  is  the  ao-reflection  region. 


Ms 


Fig.  7e  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  (M,,l,)  plane  for  vibrational  equilibrium  CO,  and 
experimental  data.  For  SMR-CMR  transition  boundnry,  solid  line 
represents  M,r*  I  and  dashed  line  represents  4*90  deg.  Dash-dot 
lines  are  for  frozen  CO,  (y  *  7/S).  Experimental  data:  o  RR,  o  SMR, 
a  CMR,  a  DMR  (Ref.  2S). 
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Fig.  7f  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  (M,,l.)  pla>te  for  CO,  ns  a  perfect  gas  (y  =  1.31) 
and  experimental  data.  For  SMR-CMR  transition  boundary,  solid 
line  represents  M,r  *  I  and  dashed  line  represents  4*90  deg.  Ex¬ 
perimental  data:  oRR.  a  SMR.  sCMR.  •  DMR  (Ref.  IS). 
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distance  between  the  two  triple  points  TT'  goes  to  zero  there. 

As  a  consequence,  the  two  DMR  runs  in  CMR  now  lie  in  the 
appropriate  regime.  Real-gas  effects  shown  by  the  equilibrium 
dashed  lines  do  not  explain  the  results,  since  the  runs  were 
done  at  different  initial  pressures.  Eight  DMR  runs  lie  in  the 
CMR  regime  and  three  SMR  runs  lie  in  the  CMR  regime.  Con¬ 
sequently,  it  appears  that  better  criteria  are  required  for  the 
rR-MR,  SMR  —CMR,  and  CMR  — DMR  transition  lines. 
Similar  remarks  apply  to  Fig.  6d.  Overall,  the  6  =  90-deg  line  is 
not  helpful  in  this  case. 

Figure  6e  includes  the  detailed  results  of  Ikui  et  al."  for  air. 

|l  is  seen  that  excellent  agreement  is  obtained  for  RR  at  the 
rR-MR  transition  line.  For  A/,  >2,  the  runs  lie  marginally 
below  it.  One  RR  run  lies  in  the  SMR  regime.  Eight  SMR’s  lie 
in  the  CMR  regime.  All  CMR  and  DMR  runs  lie  in  their  ap¬ 
propriate  regimes.  The  perfect-gas  (>  =  7/5)  transition  lines 
appear  quite  appropriate  for  their  results.  No  significant 
evidence  of  the  persistence  of  RR  into  MR  regimes  can  be 
seen.  The  dashed  line  for  6  =  90  deg  is  not  helpful  in  this  case, 
but  a  new  SMR-CMR  line  closer  to  the  CMR-DMR  line 
would  place  all  of  the  SMR  runs  in  their  appropriate  regime. 

Figures  7a-7f  deal  with  COj  in  the  (M„8m)  or  (A/,,0i) 
plane.  Figure  7a  shows  the  experiments  for  transition  lines 
based  on  a  perfect  gas,  >  =  1 .29.  Although  such  lines  are  in  er¬ 
ror,  since  the  vibrational  relaxation  lengths  are  short,  and  an 
equilibrium  flow  is  appropriate,  nevertheless,  excellent  agree¬ 
ment  is  obtained,  by  and  large,  for  all  points  except  two 
OMR’s  at  A/,-  3.5  in  the  CMR  regime  and  one  at  A/,-  4.7. 
One  CMR  at  A/,  -4  agrees  better  with  the  A/jr -  1  criterion, 
and  5  =  90  deg  is  not  helpful  in  this  case.  The  regular  reflection 
runs  at  0,  =  50  deg  are  not  a  test  of  the  RR  -  MR  transition,  as 
they  should  have  been  done  much  closer  to  or  below  that  line. 
Consequently,  bale  can  be  said  about  boundary-layer  effects 
in  this  case.  Somewhat  better  agreement  is  obtained  in  the 

(A/,,  O'*)  P|ane  sh«wn  in  7b- 
The  correct  equihbrium  transition  is  shown  in  Fig.  7c.  Here, 
the  frozen  transition  Unes  (>  =  7/5)  are  also  shown  for  com¬ 
parison.  The  agreement  is  quite  good  for  the  RR-MR  line  for 
the  frozen  case  and  poor  for  the  rest  of  the  frozen  Unes.  Eight 
CMR  runs  lie  in  the  equilibrium  DMR  regime  and  two  SMR 
runs  Ik  in  the  CMR  region.  In  this  particular  case,  the 
6  =  90-deg  line  puts  the  SMR  point  at  10  deg  in  its  proper 
place.  It  can  be  seen  from  Fig.  7a  that  the  agreement  there  is 
accidentally  better. 

Figure  7d  shows  an  (A/,,0;)  plot  of  the  experimental  runs  in 
equilibrium  C02.  Although  this  is  a  more  accurate  plot,  since 
X  does  not  have  to  be  known,  the  results  are  poorer,  since  nine 
CMR’s  find  themselves  in  the  DMR  region  and  one  SMR  oc¬ 
curs  in  the  CMR  regime.  AH  of  the  DMR  runs  are  located  cor¬ 
rectly,  as  well  as  all  RR  runs.  Again,  the  various  transition 
lines  crowd  together  near  8W  =  20  deg. 

Figure  7e  shows  the  detailed  results  of  Ikui  et  al.25  in  the 
(A/,,0.)  plane.  It  is  clear  that  there  is  no  agreement  with  COj 
treated  as  a  frozen  gas  (>  =  7/5).  However,  the  vibrational 
equilibrium  lines  show  very  good  agreement  with  the  ex¬ 
periments  over  the  entire  Mach  number  range.  Only  two  RR’s 
lie  below  the  RR-MR  line  in  the  SMR  region.  There  are  three 
DMR  runs  lying  in  the  RR  region  above  the  RR-MR  line. 
There  is  no  evidence  of  RR  persistence  of  viscous  boundary- 
layer  effects.  Four  SMR  runs  at  A/,  =  2  lie  marginally  in  the 
CMR  region;  eleven  CMR’s  lie  in  the  DMR  regime. 
Sometimes  it  is  difficult  to  distinguish  a  CMR  from  a  DMR, 
and  perhaps  this  might  account  for  some  of  the  discrepancy. 
In  this  case,  the  6  =  90-deg  line  makes  it  possible  for  several 
SMR  runs  at  8.  =  10  deg  to  lie  in  their  proper  regime. 

Figure  7f  shows  the  results  of  Ikui  et  al.”  plotted  on  the 
( A/,,0. )  plane  with  their  analytical  lines  for  a  perfect  gas 
(y  =  |.3l).  it  is  seen  that  good  agreement  is  obtained  with  all 
experimental  runs  by  and  large.  The  RR  — MR  line  does  not  fit 
the  results  as.  well  as  for  the  equilibrium  flow  (Fig.  7e)  since  14 
RR  runs  lie  in  the  DMR  region.  Only  two  RR  runs  lie  in  the 
SMR  regime  at  A/,-1.5  below  the  RR-MR  line  and  two  RR 


in  the  CMR  region.  One  borderline  DMR  run  lies  in  the  CMR 
regime  at  0*  =  20  deg.  The  Mu  =  1  line  puts  six  CMR's  in  the 
proper  regime,  unlike  6  =  90  deg;  the  preceding  results  are 
analogous  to  those  of  Ando  and  Glass51  (Fig.  7a). 

Figure  8a  shows  the  experimental  results  for  Freon-12,  from 
ikui  et  al.25  plotted  on  an  (A/,,0.)  plane,  for  their  transition 
lines  based  on  a  perfect  gas  (»-=  1. 141).  Again  the  perfect-gas 
lines  give  reasonable  agreement.  Four  RR  runs  lie  in  the  CMR 
region.  Many  DMR  runs  lie  in  the  CMR  regime.  About  four 
SMR’s  lie  marginally  in  the  CMR  region  and  two  CMR  runs 
lie  in  the  DMR  regime.  The  dashed  line  is  an  experimental  best 
fit  for  the  RR-TDMR.  Unfortunately,  Ikui  et  al.  did  not  ex¬ 
tend  the  TDMR  runs  to  higher  A/,  and  lower  0.  to  provide  a 
better  picture  of  the  TDMR  regime,  which  can  be  predicted 
analytically  (see  Fig.  8b).  Only  in  the  CMR  region  is  there 
evidence  of  the  persistence  of  RR.  An  equilibrium  (A/,,0.) 


Fig.  Si  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  In  the  (A/,,0.)  plane  for  Freon-12  as  a  perfect  gas  with 
>=1.141  and  experimental  data.  Dashed  line  represents  the  ex¬ 
perimental  RR-TDMR  transition  boundary.  For  SMR-CMR,  the 
solid  line  represents  Mjy  =  l,  and  the  dash-dot  line  represents  6  =  90 
deg.  Experimental  data:  TDMR  (terminal  double-Mach  reflection). 
oRR,  oSMR,  *  CMR,  a  DMR  (Ref.  25). 


Fig.  8b  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  (A/,,0.)  plane  for  Freon-12  and  experimental  data. 
Dash-dot  lines  represent  the  frozen  flow  (>  =  4/3).  Solid  lines  repre¬ 
sent  equilibrium  Freon-12.  Dash-double  dot  line  represents  the 
RR-TDMR  experimental  limit  a  TDMR  (terminal  double-Mach 
reflection).  oRR,  a  SMR,  a  CMR,  a  DMR  (Ref.  25). 
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plot  is  shown  in  Fig.  8b.  It  also  includes  the  transition  lines  for 
a  frozen  flow  with  7  =  4/3,  for  which  the  agreement  with  ex¬ 
periment  is  poor.  This  time,  more  DMR  runs  lie  in  the  RR 
region  including  the  TDMR  runs.  Note  that  the  condition 
x'  =0  yields  only  a  small  TDMR  region  below  the  RR  -  MR 
line.  Consequently,  one  can  only  suspect  that  the  detachment 
criterion  of  62  =  6^,  due  to  von  Neumann  is  invalid  for 
pseudostationary  flow.  The  equilibrium  SMR-CMR  is  not 
improved  from  the  one  in  Fig.  8a,  nor  does  the  criterion  6  =  90 
deg  help.  The  CMR-DMR  line  does  not  improve  overall  the 
location  of  the  runs.  Some  CMR  runs  are  more  poorly  located 
and  some  DMR  runs  are  better  located  in  their  respective 
regions. 

Figures  9a  and  9b  show  the  results  for  SF6  in  the  (A/,,6, ) 
and  (A/,,6^)  planes. 58  Considerable  effort  was  made  in  this 
case  to  check  the  transition  lines  themselves.  Figure  9a  shows 
the  results  for  frozen  (7  =  4/3)  and  equilibrium  flows.  The 
frozen-flow  lines  do  not  represent  the  experiments.  Two  RR 
and  two  DMR  runs  lie  in  the  CMR  region  below  the  RR  — MR 
transition  lines.  Six  RR  runs  lie  right  on  the  RR  — MR  transi¬ 
tion  line.  There  is  no  serious  evidence  of  RR  persistence  up  to 
A/,-6.5  and  two  TDMR’s  at  A/,  -  8.  One  CMR  run  lies  in  the 


SMR  region  and  five  DMR’s  lie  in  the  CMR  regime.  The 
A/2r  =  1,  rather  than  6  =  90  deg,  puts  five  CMR  runs  in  their 
proper  place.  However,  overall  the  agreement  is  good.  Similar 
remarks  can  be  made  for  the  (A/,,6^)  plot  of  Fig.  9b.  Here, 
however,  it  is  clear  that  DMR  persists  into  the  RR  regime  over 
the  entire  range  of  A/,.  Consequently,  the  (A/,,6,)  plot  gives 
poorer  rather  than  better  agreement  with  the  RR  — MR  transi¬ 
tion  line.  It  should  be  noted  that,  if  three-shock  theory  was 
used  here  rather  than  two,  x  would  have  a  value  from  1  to  3 
deg  over  the  RR  — MR  line,  which  would  be  pushed  upward 
for  better  agreement  with  experiments  than  the  present  line. 
The  small  dashed  line  running  from  CMR-DMR  to  point  P 
on  the  RR— MR  transition  line  indicates  the  fact  that  in  ac¬ 
tuality  this  takes  place.  It  makes  for  better  agreement  in  Figs. 
9a  and  9b  in  that  the  DMR  runs  at  A/,  -  2  lie  in  their  correct 
region.  Again,  the  transition  lines  are  too  close,  especially  at 
6,  -  10  deg  in  Fig.  9b. 

Figure  9c  shows  the  SF6  results  in  the  (A/,,6,)  plane  for  a 
perfect  gas  (7=1 .093).  The  experiments  are  quite  well 
represented  on  this  plot.  Two  RR  runs  lie  below  the  RR-MR 


Fig.  9a  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  (A/,,!,)  plane  for  vibrational  equilibrium  SF,  and 
experimental  data.  For  SMR-CMR  transition  boundary,  solid  line 
represents  M]r  =  l  and  dashed  line  represents  6  =  90  deg.  Dash-dot 
lines  are  for  frozen  SFt  (7  =  4/3).  Experimental  data:  a  TDMR, 
oRR,  o  SMR,  a  CMR.  •  DMR  (Ref.  3*). 


Fig.  9c  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  (A/,,6,)  plane  for  SF,  as  a  perfect  gas  (7=  1.093) 
and  experimental  data.  For  SMR-CMR  transition  boundary,  solid 
line  represents  M!r  =  1  and  dashed  line  represents  6  =  90  deg.  Ex¬ 
perimental  data:  aTDMR,  oRR,  oSMR,  a  CMR,  •  DMR  (Ref. 
38). 


Fig.  9b  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  )  plane  for  vibrational  equilibrium  SFt  and 

experimental  data.  For  SMR-CMR  transition  boundary,  solid  line 
represents  M]r=  I  and  dashed  line  represents  6  =  90  deg.  Dash-dot 
lines  are  for  frozen  51,  (7  =  4/3).  Experimental  data:  aTDMR, 
oRR.  oSMR,  -a CMR,  •  DMR  (Ref.  38).  NR  Is  the  no-reflection 
region. 


Fig.  9d  Regions  and  transition  lines  of  regular  (RR)  and  Mach  (MR) 
reflections  in  the  (A/,,6, )  plane  for  SF6  as  a  perfect  gas  (>  =  1.093) 
and  experimental  data.  For  SMR-CMR  transition  boundary,  solid 
line  represents  M2r=l  and  dashed  line  represents  6  =  90  deg.  Ex¬ 
perimental  data:  aTDMR.  oRR,  oSMR,  a  CMR.  •  DMR  (Ref. 
38).  NR  is  the  no-reflection  region. 
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Fig.  10a  TDMR  regions  la  ( M,,tw )  plane  for  vibrational  equilibrium 
Freon-12  and  lor  perfect  gases  with  y  =  1.10,  1.093,  and  1.08. 


Ms 


Fig.  10b  TDMR  regions  in  (A #,,#„)  plane  for  vibrational  equilibrium 
SFt,  f-1.06,  and  >- 1.04. 


2  4  M  6  8  10 


Fig.  11  Variation  of  departure  of  experimental  5,from  frozen-gas 
calculation  (J,  —  hj)  normalized  by  the  difference  between 
equilibrium-  and  frozen-gas  calculations  (4,  -hj)  (Ref.  43). 


transition  line  in  the  CMR  region  and  four  below  in  the  DMR 
regime  at  A#, -2,  and  four  below  at  A/,  -4,  showing  a  per¬ 
sistence  of  RR.  Several  DMR  runs  lie  in  the  CMR  region  at 
0„~  45,  43,  40,  20  deg,  and  especially  at  10  deg.  All  CMR 
runs  lie  properly  in  their  region  at  or  below  the  M1T-  1  line, 
except  two  in  the  SMR  region  at  A/,-3.5.  The  equilibrium 
plot  of  Fig.  9a  is  superior  since  the  RR  — MR  line  lies  lower  for 
lower  y  than  for  the  fixed  (7  =  1.093)  line  of  Fig.  9c.  The 
SMR  — CMR  and  CMR  — DMR  lines  are  also  better  in  Fig.  9a 
for  the  same  reason,  providing  better  agreement  with  experi¬ 
ment  overall  for  the  equilibrium  plot. 

Figure  9d  provides  the  (A/,,0;)  plot  for  7=  1.093.  For 
overall  agreement,  it  is  not  as  good  as  Fig.  9b.  Here  three 
DMR's  lie  "in  the  CMR  region  at  A/,  -  3.5  and  several  at 
A/,-2.  These  would  lie  correctly  in  the  DMR  region  if  the 
CMR— DMR  line  were  curved  to  join  the  SMR  — CMR  line  to 


2  4  6  8  10 

M, 


Fig.  12  Variatiou  of  experimental «  with  incident  shock-wave  Macb 
number  M,  (Ref.  43):  a ,  N,  (Ref.  30);  ■ ,  Ar  (Ref.  30);  • ,  CO, 
(Ref.  37);  v  ,  air  (Ref.  36). 
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Fig.  13  Difference  In  the  incident  wave  angle  to  the  Mach  stem  44, 
between  experimental  and  numerical  results  for  CO,  (Ref.  43):  Ex¬ 
perimental  mails  from  Ref.  37. 


the  RR  —  MR  line  at  point  P,  as  noted  previously.  On  the  other 
hand,  several  DMR’s  and  the  two  TDMR’s  lie  in  the  RR 
region.  Overall,  the  equilibrium  curves  of  Fig.  9b  provide  bet¬ 
ter  agreement. 

Figure  10a  shows  the  extent  of  the  TDMR  over  a  range  of  7 
as  well  as  for  vibrational  equilibrium  Freon- 12  and  SFt.  It  is 
seen  that  the  regions  are  bounded  by  the  two-shock  theory 
RR-MR  transition  line  and  the  line  along  which  \  '  =0.  The 
region  for  a  perfect  gas  with.y  =  1. 10  is  very  small,  and  occurs 
for  A/,  >9  and  0„  -  38  deg.  For  a  perfect  gas  with  7=  1.093, 
the  region  is  much  larger  and  occurs  for  A/,  >7.5  and  0„>37 
deg.  For  vibrational  equilibrium  Freon- 12,  the  region  is  much 
larger  still  and  starts  at  A/,  -  7.  For  7  =  1 .08,  the  TDMR  region 
is  quite  large  and  starts  at  A/,  -  5.  Note  that,  unlike  the  results 
of  Ikui  et  al.  shown  in  Fig.  8a,  the  TDMR  experiments  are 
supposed  to  lie  in  the  DMR  region,  not  in  the  RR  regime,  as 
they  found.  It  is  not  clear  why  this  has  happened.  Even  if  the 
three-shock  theory  RR  —  MR  line  was  used  and  x  ’  was  in  error 
by  1  or  2  deg,  it  would  not  account  for  this  larger  discrepancy, 
as  can  be  seen  from  Fig.  10b.  One  can  only  speculate  that  the 
RR-MR  transition  line  is  not  appropriate  for  pseudosta¬ 
tionary  flows. 

Figure  10b  shows  the  regions  of  TDMR  for  vibrational 
equilibrium  SFt,  7=  1.06  and  1.04.  It  is  seen  that  the  regions 
become  progressively  larger.  It  should  be  possible  to  observe 
TDMR  in  vibrational  equilibrium  SF6  at  A/, -5.  Ikui  et  al.25 
observed  TDMR  in  Freon- 12  at  A/,  -  3.8  and  0„  -46  deg,  and 
the  present  analysis  shows  this  to  occur  at  M,  -  7  and  0.  -  37 
deg.  Consequently,  this  problem  will  require  a  resolution  in 
the  near  future. 

In  summary,  it  can  be  stated  that  the  (A/„0„)  plots  provide 
the  experimenter  with  good  engineering  results  for  all  gases. 
Whether  a  point  lies  in  one  region  or  another  sometimes 
depends  on  personal  judgment  when  the  result  is  not  dear  cut. 
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Consequently,  it  is  more  important  to  make  measurements  of 
the  angles  4,  and  u' ,  which  are  quantitative  experimental 
results,  and  compare  them  with  their  counterparts  from 
analysis  using  two-  and  three-shock  theories,  as  applicable,  in 
order  to  provide  better  tests  for  the  state  of  the  gas— frozen  or 
equilibrium.  The  induced  boundary  layer  causes  5  to  be  larger 
and  u’  to  be  smaller.  Consequently,  one  can  learn  a  good  deal 
about  the  boundary  layer  itself  from  such  measurements.  The 
angle  x  can  be  affected  by  the  fact  that  the  inclination  of  the 
slipstream  S  depends  on  the  two  velocities  on  either  side  of  the 
ideal  contact  surface,  one  in  state  (2)  moving  at  a  higher 
velocity  than  that  in  state  (3).  Consequently,  the  velocities  in 
states  (2)  and  (3)  diverge  due  to  the  layer  thickness,  giving  rise 
to  a  deflection,  a,  which  could  affect  x.  as  well  as  the  relation 
4,  -42  =  4Jf  used  in  the  three-shock  theory. 

Critique 

It  has  been  shown  above  that  many  experiments  in  various 
gases  have  now  been  performed  on  regular  and  Mach  reflec¬ 
tions  of  oblique  shock  waves  in  pseudostationary  flow.  The 
experimental  agreement  with  the  analytical  boundaries  for 
such  reflections  using  two-  and  three-shock  theories  are  rea¬ 
sonable  but  not  precise  enough  over  the  entire  range  of  inci¬ 
dent  shock-wave  Mach  numbers  (M,)  and  compression  wedge 
angle  (<?„).  In  order  to  improve  the  agreement,  the  assump¬ 
tions  and  criteria  employed  in  the  analysis  were  critically  ex¬ 
amined  using  the  foregoing  experimental  data.4’ 

Several  criteria  were  proposed  to  predict  analytically  the 
transition  boundaries  between  the  various  reflections. 
However,  the  first  complete  transition  solution  in  the  (ms,8w) 
or  (M„8’w)  plane  using  two-  and  three-shock  theories  and  the 
relevant  criteria  to  obtain  the  transition  lines  were  given  in 
Refs.  30  and  31,  respectively,  for  N2  and  Ar;  in  Ref.  37  for 
C02;  in  Ref.  28  for  perfect  and  equilibrium  air  up  to  M,  =  20; 
in  Ref.  38  for  SF4;  and  in  the  present  paper  for  Freon-12. 23 

Comparison  of  Experimental  and  Numerical  Values  of  4 

It  has  been  noted  previously  that  by  measuring  the  angles  6, 
X,  and  b>',  it  is  possible  to  obtain  quantitative  checks  with 
analysis  as  well  as  inferring  the  state  of  the  gas,  whether  it  was 
frozen  or  in  equilibrium.  A  minimum  finite  length  of  1  mm  is 
required  to  measure  the  slope  of  a  shock  wave  from  an  in- 
terferogram.  The  measured  angle  is  then  the  average  slope 
within  this  length.  The  vibrational  relaxation  length  tv  at 
P0  -  15  Torr  and  T0  =  300  K  becomes”  I  mm  at  M,  -  8  for  02, 
M,>  10  for  N2,  A/,-5  for  C02,  and  M,~  2  for  SF4.  In  two 
and  three  shock-wave  systems  for  regular  and  Mach  reflec¬ 
tions,  the  relaxation  lengths  do  not  vary  significantly43  com¬ 


mit.  14  Compirisoa  of  experimental  reflected  wave  angle  w  with 
numerical  results  for  various  displacement  angles  for  C02:  Ex¬ 
perimental  points  (Ref.  37):  v -45.5  deg:  a,  0„-47  deg; 
o,#„-50  deg;  4  -  40  deg;  lines  (equilibrium-gas  calculation!: 

- 4,-0  deg; - .  -  I  deg;  —  -  2  deg  (Ref.  431. 


pared  to  a  I0M07  variation  in  f„  for  1<M,<\0,  Conse¬ 
quently,  the  Mach  number  at  which  f„  =  I  mm  behind  each 
shock  wave  in  Mach  reflection  does  not  differ  very  much. 
Therefore,  the  choice  of  1  mm  as  a  characteristic  flow  length  is 
quite  reasonable.  When  the  relaxation  length  is  of  the  same 
order  as  the  characteristic  length,  the  simplified  two-  and 
three-shock  theories  can  no  longer  be  used.  Therefore,  it  is 
reasonable  to  apply  a  frozen  ((„►!  mm)  flow  or  an 
equilibrium  (f„  <  1  mm)  fle  w  assumption.  Only  vibrational  ex¬ 
citation  need  be  considered  in  the  present  experiments  as 
dissociation,  electronic  excitation,  and  ionization  relaxation 
lengths  are  much  longer  than  1  mm  in  the  present  shock  Mach 
number  and  pressure  ranges  of  the  experiments.  Air  with  21  Oh 
02  is  an  exception  at  high  shock  Mach  numbers.  However, 
dissociation  is  negligible  for  M,  <  10.  The  specific  heat  ratio 
for  frozen-gas  calculation  is  5/3  for  Ar;  7/5  for  N2,  02,  air, 
and  C02  (linear  molecule);  and  4/3  for  Freon- 12  and  SF4. 

In  Ref.  43,  there  are  several  plots  of  the  variation  of  angle  4, 
with  shock  Mach  number  for  fixed  effective  wedge  angle  B'm 
for  C02,  N2,  and  Ar,  and  in  Ref.  38  for  SF4.  These  may  be 
condensed  by  plotting  the  quantity  *  =  (4, -6/)/(6, -6f) 
where  4,  is  the  measured  4,  6f  the  calculated  frozen  4,  and  4, 
the  calculated  equilibrium  4.  It  is  seen  from  Fig.  11  for  C02, 
for  example,  that  for  M,  >4,  k- 1;  i.e.,  4,  =4,  and  the  agree¬ 
ment  is  for  C02  in  equilibrium,  despite  the  fact  that  the  ex¬ 
periments  in  the  ( M,.8W )  plane  agree  best  with  a  perfect  gas 
(7=  1.29).  These  results  are  more  reliable:  As  no  transition 
criteria  are  required  in  the  calculation  of  4,  a  quantitative  com¬ 
parison  can  be  made  for  each  pair  of  analytical  and  ex¬ 
perimental  points;  a  measurement  of  4  has  no  ambiguity 
unlike  the  classification  of  the  reflection  types  near  the  transi¬ 
tion  boundaries  where,  for  example,  a  CMR  may  look  like  an 
SMR  or  a  CMR  like  a  DMR.  The  parameter  4^  should  be  used 
instead  of  Bm  in  the  comparisons,  since  8 'w  can  be  obtained 
directly  from  Eqs.  (1-4),  while  8W  is  derived  by  assuming  that 
the  Mach  stem  M  is  perpendicular  to  the  wedge  surface  at  the 
triple  point.  This  is  not  always  the  case  experimentally,  as 
shown  subsequently. 

Effects  of  Mach  Stem  Curvature 

In  Figs.  1  and  2,  all  of  the  Mach  stems  M  are  drawn  ideally 
perpendicular  to  the  wedge  surface  from  the  triple  point  T.  In 
reality,  the  Mach  stem  can  have  a  convex  curvature  ( + )  at  an 
angle  1  from  T  or  with  a  concave  ( -  )  curvature  but  again  end¬ 
ing  up  perpendicular  to  the  wedge  surface. 


a!  Intcrferogram  and  experimental  hopycaks  region  and  (0) 
1 .00;  (I)  2.71 ;  (2)  3.68;  (3)  3.33;  (a)  3.36;  (b)  3.44;  (c)  3.52;  (d)  3.60;  (c) 
3.68;  (0  3.75;  (g)  3.83;  (b)  3.91;  (i)  3.99;  <J)  4.06;  (k)  4.14. 


b)  Numerically  simulated  isopycniev  by  Glu  el  al.14  shown  in  f'ig. 
15a. 


f  ig.  15  SMR.  \f,  =  2.03,  #.  =  27  deg.  p,  -  33.3  hPa.  f,  =  299.2  K. 
Air,  ?4  =  0.387  kg  m3. 
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The  Mach  stem  was  assumed  perpendicular  to  the  wedge 
surface  at  the  triple  point  in  deriving  6,  from  4,  in  the  three- 
shock  theory.  The  deviation  t  from  this  idealization  is  shown 
for  Ar.  Nj,  air,  and  C02  in  Fig.  12.  It  is  seen  that  t  has 
positive  values  (0-8  deg)  for  M,  >4  and  decreases  to  negative 
values  ( -  3  deg)  at  M ,  =  2.  Consequently,  6W  as  derived  from 
0,  has  an  error  equal  to  t.  Therefore,  comparisons  of  quan¬ 
tities  in  Mach  reflection  should  be  made  with  B'w,  as  long  as  i 
cannot  be  predicted  accurately.  This  error  is  the  main  reason 
for  the  inaccurate  prediction  of  *  using  the  three-shock 
theory. 

In  Fig.  13,  the  effect  of  <  on  4,  is  shown  for  C02  as  the 
variation  of  44,  =  4 ,  (numerical)  -  4,  (experimental)  with 
shock  Mach  number  M,  The  numerical  results  are  based  on 
an  equilibrium  gas,  verified  in  Fig.  II.  It  is  seen  that  44,  has  a 
value  of  about  -  2  deg,  independent  of  M,.  The  results  for 
other  gases  are  similar  albeit  somewhat  more  dispersed. 

llffd  of  Slipstream  Thickness 

In  the  solution  of  Eqs.  (1-4),  only  the  pressures  (Eq.  (7)J 
and  flow  directions  [Eq.  (8)]  were  assumed  to  be  continuous 
across  the  slipstream,  an  idealized  surface.  In  reality,  it  is  a 
shear  layer  of  increasing  finite  thickness,~where  the  velocity  in 
state  (2)  is  greater  than  that  in  state  (3)  and  the  temperature  in 
state  (3)  is  greater  than  that  in  state  (2).  Near  the  triple  point 
the  layers  are  laminar;  beyond  they  become  turbulent.  Conse¬ 
quently,  as  noted  previously,  the  shear-layer  thickness  gives 
rise  to  a  deflection  of  the  velocities  in  states  (2)  and  (3).  If  it  is 
assumed  that  the  velocity  displacement  in  the  layer  a  exists  at 
the  triple  point,  then  its  possible  affect  on  x  should  be  taken 
into  account  in  the  discrepancies  between  the  numerical  and 
experimental  values  of  x- 

Effect  of  Shock-Induced  Boundary-Layer 
on  the  Wedge  Surface 

It  has  been  found  by  many  experimenters  that  RR  persists 
into  the  MR  regime.  This  has  been  called  the  "von  Neumann 
paradox."  Recently,  several  analyses19,4345  have  shown  that 
this  is  due  to  the  boundary-layer  negative  displacement 
thickness.  Consequently,  the  displaced  wall  is  below  the  actual 
wall.  The  reflected  wave  R  moves  below  to  a  new  position  so 
that  the  angle  6  between  1  and  R  is  now  larger  or  the  reflected 
wave  angle  u>'  is  smaller.55  The  condition  8,  -82=8d  must 
now  be  used  instead  of  Eq.  (6).  This  can  be  checked  and  has 
been  verified  experimentally.  Figure  14  compares  experimen¬ 
tal  w'  for  C02  with  calculations  for  several  values  of  8d.  The 
experimental  results  lie  between  the  lines  of8d-0  and  -  2  deg. 
Similar  results  were  obtained  for  Ar,  N2,  and  air.  A  displace¬ 
ment  angle  8d  =  -  I  deg  shifts  the  RR-MR  von  Neumann 
detachment  transition  boundary  by  40.=  -0.7  deg  for  N2 
and  air,  and  40.  =  -0.5  for  C02.  Recently,  Sakurai44  pro¬ 
posed  to  use  a  pseudostationary  RR-MR  criterion.  However, 
its  effectiveness  is  yet  to  be  tested  experimentally. 

Flowfield  Solutions 

The  foregoing  dealt  with  simplified  algebraic  solutions**  for 
RR  and  MR,  which  involved  mainly  the  shock  waves  and  the 
slipstream.  In  order  to  compare  experimental  and  analytical 
solutions  of  the  entire  flowfield  for  such  sensitive  parameters 
as  the  density  distribution,  which  can  be  measured  inter- 
ferometrically,  it  is  necessary  to  solve  the  time-dependent”  M 
or  pseudostationary  Euler  equations  as  noted  subsequently, 
for  a  perfect  or  equilibrium  gas.  Otherwise,  it  would  be 
necessary  to  use  rate  equations  for  equilibration  as  well. 

Continuity: 

p,  +  (pit),  +  (pc),  -  0 

Momentum: 

(pu),  +  {pul  +p),  +  (pur),  =0 

(pt  ),  +  (pi  ),  +  (pi-’  +  />)„  -  0 


Energy:  j 

{pE),+ (puE +  up),  +  (pvE+ vp),=0  i 

i 

where  p  is  the  density,  u=(u,t>)  is  the  velocity  field,  ! 

E  =  ‘/j(u2  +  v1)  +  e  is  the  total  specific  energy,  e  is  the  specific  1 

internal  energy,  and  p  is  the  pressure.  The  system  is  closed  by  I 

specifying  an  equation  of  state  (EOS).  ; 

« 

p  =  p(p.t)  (10)  ] 

Such  an  equation  can  be  obtained  from  the  thermal  equation 

of  state  p=pRT  and  the  specific  internal  energy  for  a  perfect  ] 

gas  e  =  C,T,  to  yield  j 

p=(y-l)pe  (II) 

where  y  >  1  is  the  ratio  of  specific  heats  and  is  considered  to  be 
a  constant.  For  gases  such  as  C02  and  SFt  in  vibrational 
equilibrium  the  equation  of  state  has  to  be  modified.  For  self-  I 

similar  motion,  the  problem  has  no  intrinsic  length  scale.  Con-  | 

sequently,  a  new  pseudostationary  coordinate  system  can  be 
used,  such  that  ({,>»)  =  I  (jt — jr0 ) /(/ - /0).  0'->o)/('-'o)). 
where  (x0.  >o)  are  the  coordinates  of  the  wedge  corner  C,  and 
t0  is  the  time  when  the  incident  shock  wave  /  reaches  the  cor¬ 
ner.  As  shown  by  Jones  et  al.,21  the  system  of  Eq.  (9)  can  be 
transformed  into  a  pseudostationary  system  [Eq.  (12)]  .M 

(pw),  +  (ptf),=  -2p  | 

{pu2  +p){  +  (put>),=  -3pu 
(puv)l  +  (pi? +/>),=  -ipv 

I 

(pu//){  +  (pt5/f),=  -p(w2  +  t?)-2pH  (12)  I 

where 

«  =  «-{,  =  H=  H(£2  +  (?)  +  *  (13) 

and  h=e  +  p/p  is  the  specific  enthalpy.  Also  ( u,v )  and  //are 
referred  to  as  the  self-similar  velocity  field  and  self-similar 
total  enthalpy,  respectively.  In  addition  M  is  defined  by 

\f  =  (iJ  +  i?)/a2  (14) 

where  a  is  the  sound  speed  and  Af  is  called  the  self-similar 
Mach  number.  The  system  of  Eq.  (12)  is  the  steady  Euler 
equations  with  the  addition  of  source  terms.  Note  that  the 
ratio  s/L  is  constant  (Fig.  Ic)  for  given  initial  conditions,  for 
self-similar  solutions  of  the  nonstationary  equations  just  as  j  is 
contant  for  steady  supersonic  flow.  In  this  and  other  ways,  a 
change  to  pseudostationary  coordinates  is  very  useful  in  the 
analysis  of  such  flow  fields  and  was  used  in  numerical  simula¬ 
tions  of  the  interferometric  experiments  on  oblique  shock- 
wave  reflections.'’ 

Perhap1-  the  best  example  of  showing  the  great  strides  that 
have  taken  place  in  the  reliability  of  numerical  simulation  of 
such  problems  is  to  compare  the  results  of  Schneyer,*'  of  a 
decade  ago,  with  the  present  state-of-the-art  simulation  by 
Glaz  et  al.*4  Figure  15a  shows  an  interferogram  of  an  SMR  in 
air  at  M,  =  2.03,  8m  =  27  deg,  p0  =  250  Torr,  and  T0  =  300  K. 

The  lines  of  constant  density  (isopycnics)  over  the  flowfield 
are  also  shown  in  the  caption.  Figure  15b  shows  the  same  case 
treated  by  Glaz  et  al.’4  a  decade  later,  which  is  a  very  good 
simulation  indeed.  Many  such  cases  were  treated  in  air,  Ar, 
and  SF*.’4"  49  The  agreements  with  experiments  are  all  sur¬ 
prisingly  good  considering  that  the  simulations  utilized  the 
l  uler  equations  for  inviscid  flow  and  the  inierferogtams  are 
for  real  flows  with  viscosity  and  real-gas  effects.  Needless  to 
say,  the  simulations  provide  numerous  data  on  pressure,  tem¬ 
perature,  Mach  number,  enthalpy,  etc.,  that  cannot  readily  be 
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Fig.  16a  Collision  of  OMR  with  a  90-deg  ramp  la  CO, ,  I.  -  20 deg,  f -  134  nun,  A  - 12.7  ana,  /*,-!•  Ton,  A - 6943  A.  Aa  interferograia  taken 
al  /  -  56  a*.  M,  -  Ml,  7*  -267  K,  aad  e0  ”  2-3*  x  10  '  r  g/c«,J.  A*/>»  -  1.27. 1,  incident  ihock  wave;  M,  Dm  Mack  stem;  M' ,  second  Mack  stem; 
ft,  llnl  reflected  skock  wave;  ft  ,  second  reflected  skock  wave;  T,  first  triple  point;  T  second  triple  point  (Ref.  24). 


Fig.  16k  Aa  Inlerferogram  taken  ill*  1 19  ns,  M ,  «  5.78,  Tt  »  294  K,  ■  2.40  x  10  ’  g  cm1,  Wi  =  1-26.  B.  bif  area  lion  of  shock  wave;  C(,  con¬ 
tact  front;  £,  espansioa  wave;  /,  Incident  shock  wave;  P,  protnistion  doe  to  local  sidewall  boundary -layer  separation;  .V,to  .V?,  shock  waves;  ft, 
rarefaction  wave. 
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Fig.  17a  Numerical  simulation  of  the  collision  of  a  DMR  with  a 
00 -de(t  ramp  in  C02-  Thirl)  equally  spaced  numerical  isopycnics. 
Conditions  as  for  Tig.  16a.  Note  that  the  interaction  of  the  Mach  stem 
with  the  rolled-up  slipstream  is  more  severe  than  in  the  interferogram. 
In  addition,  a  second  triple-point  system  and  slipstream  appear  at  the 
kink  in  the  Mach  stem.  A  second  normal  shock  wave,  to  the  left  of  the 
curved  Mach  stem,  separated  by  a  contact  surface,  can  also  be  seen. 
These  transitions  do  not  appear  in  the  interferogram  and  are  sup¬ 
pressed  by  the  induced  boundary  layer.  (See  Refs.  33  and  34  for 
similar  results  in  air,  etc.). 


Figs.  17b  Numerical  simulation  of  the  collision  of  a  DMR  with  a 
90-deg  ramp  in  COt.  Thirty  equally  spaced  numerical  isopycnics. 
Conditions  as  for  Fig.  16b.  The  inviscid  results  are  modeled  quite  well 
except  for  the  position  of  contact  surface  C(  and  the  recompression 
shock  wave  ij.  However,  the  viscous  effects,  such  as  the  shock-wave 
bifurcation  B  and  the  protrusion  /*,  due  to  the  local  boundary-layer 
separation  on  the  window,  are  not  modeled,  as  expected.  This  effec¬ 
tively  points  out  the  necessity  of  redoing  such  simulations  using  the 
Navier-Slokes  equations. 


measured  experimentally.  Undoubtedly,  repeat  simulations 
using  the  Navier-Stokes  equations  and  equilibration  rate  equa¬ 
tions  for  the  degrees  of  excitation  involved  would  provide 
some  important  answers  as  to  whether  the  costly  refinements 
were  warranted. 

It  is  worth  noting  that  some  very  worthwhile  concepts  have 
been  added  to  pseudostationary  oblique  shock-wave  reflec¬ 
tions  recently,  using  the  Monte-Carlo  technique.  Here  the 
reflections  are  simulated  at  the  molecular  level  by  using  a  large 
number  of  particles  ( >  10,000)  and  following  their  positions  in 
phase  space.50'51  Of  the  two  papers,  Seiler’s  paper51  is  of 
greater  interest  since  he  treated  both  specular  and  diffuse 
reflections  of  the  hard-sphere  molecules  and,  therefore,  was 
able  to  obtain  the  boundary  layer  induced  by  the  incident  and 
reflected  shock  waves  in  RR  at  the  reflection  point.  Seiler’s 
results  confirm  that  RR  persists  into  the  MR  regimes  as  a 
result  of  the  induced  boundary  layer.  He  also  showed  that  the 
leflected  wave  angle  w'  is  smaller  as  a  result.  However,  the 
method  will  have  to  be  refined  to  produce  sharper  shock 
waves,  well-defined  slipstreams,  and  more  meaningful  isolines 
compared  with  experimental  interferograms  and  the  numerical 
simulations  obtained  by  Glaz  et  al.n  ,5-49  Unfortunately,  Seiler 
did  not  investigate  the  expansion-compression  wave  system 
model  produced  by  the  boundary-layer  displacement  thickness 


at  its  leading  edge,  used  by  several  investigators  in  order  to 
substantiate  if  at  the  molecular  level  the  model  is  correct,  since 
it  cannot  be  seen  in  optical  enlargements  of  the  flow. 

A  more  complex  simulation  is  that  of  a  problem  where  self¬ 
similarity  cannot  be  applied.  For  example,  the  collision  of  a 
Mach  reflection  with  a  90-deg  ramp  on  the  wedge  surface  pro¬ 
vides  such  a  case.24  4*  Figure  16a  shows  an  interferogram  of  a 
DMR  in  COj .  Ms  =  5.81 . 6W  =  20  deg,  p0  =  10  Torr,  T0  =  297  K 
before  colliding  with  a  12.7-mm  ramp.  Figure  16b  shows  the 
wave  system  60  its  later,  after  the  collision  had  taken  place.  It 
is  seen  that  the  wave  system  involving  Mach  stem-slipstream 
interactions,  boundary-layer  interactions,  rarefaction,  and  ex¬ 
pansion  waves,  etc.,  has  become  very  complex.  The  type  of  in¬ 
teraction  is  dependent  on  the  relative  heights  of  Mach  stem 
and  ramp.  Figure  17  shows  the  numerical  simulations.  Sur¬ 
prisingly  good  agreement  is  obtained  of  the  main  features  of 
the  flow  using  the  Euler  equations  of  motion.  However,  the 
viscous  effects  are  not  simulated  and  effectively  show  the  need 
for  using  the  Navier-Stokes  equations  for  adequate  simula¬ 
tion.  Many  additional  details  can  be  found  in  Refs.  24  and  49. 

Conclusions 

1)  Although  much  has  been  learned  about  oblique  shock- 
wave  reflections  since  Mach22  wrote  his  pioneering  paper  in 
1878,  a  great  deal  is  still  to  be  known.  The  best  available  data 
to  date  from  several  sources  on  such  reflections  have  been 
presented  graphically  in  an  consistent  format. 

2)  The  presently  accepted  transition  lines  for  RR -MR, 
SMR-CMR,  and  CMR— DMR  provide  very  reasonable 
engineering  answers  as  to  the  type  of  reflection  that  will  occur, 
given  M,,6W ,  and  y  for  frozen  and  perfect  gases  and  for  vibra¬ 
tional  equilibrium  given  additionally  T0.  For  gases  with  low  y 
at  room  temperature,  such  as  COj,  Freon-12,  and  SF6,  real- 
gas  effects  in  vibration  are  very  important  and  must  be  taken 
into  account.  They  are  not  precise  enough,  and  improved 
transition-line  criteria  will  have  to  be  found  to  account  for 
those  experiments  that  give  rise  to  the  "von  Neumann 
paradox”  for  the  RR-MR  transition.  For  the  SMR-CMR 
and  CMR -DMR  lines  at  lower  dw,  where  a  number  of  ex¬ 
periments  fall  outside  the  present  lines  (see  Figs.  7-9),  new 
criteria  will  have  to  remedy  this  situation. 

3)  The  experimental  CMR  -  DMR  line  leading  to  point  P  on 
the  RR  -  MR  line  will  have  to  be  replaced  by  an  analytical  line 
(see  Figs.  9a,  9b,  and  9d). 

4)  It  is  not  clear  whether  the  induced  boundary-layer  slope 
or  thickness  is  the  important  parameter  to  account  for  the 
"von  Neumann  paradox."  Some  experiments  with  induced 
turbulent  boundary  layers  caused  by  specified  roughness 
should  help  in  a  decision  on  this  point. 

5)  It  is  not  known  why  the  TDMR  experiments  lie  in  the  RR 
regime  rather  than  in  the  DMR  region  as  predicted  from  the 
condition  x'  =0-  Perhaps  the  von  Neumann  RR  — MR  line  is 
not  applicable  to  pseudostationary  flow. 

6)  More  needs  to  be  known  about  RR  and  MR  over  rough 
surfaces,  spongy  surfaces,  surfaces  with  a  thin  layer  of  He  or 
other  gases,  and  in  a  dusty-air  environment. 

7)  Improved  computational  simulations  using  the  Navier- 
Stokes  equations  and  real-gas  rate  equations  are  required  for  a 
comparison  with  the  available  interferometric  data. 

It  can  be  seen  that  much  work  still  remains  to  be  done  in 
pscudostationary  oblique  shock-wave  reflections  in  order  to 
answer  some  of  the  above-noted  problems.  Undoubtedly,  such 
work  will  be  forthcoming  in  the  near  future.  In  addition  to  the 
single-wedge  problems,  the  degree  of  complexity  increases 
using  multiple  straight  wedges,  as  well  as  convex  and  concave 
wedges. 5)'54  For  further  details  see  Ref.  55. 
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